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NEW GENERALIZED RHEOLOGICAL MODEL FOR LUBRICATION OF A 
BALL SPINNING IN A NONCONFORMING GROOVE 
by Charles W. Allen, Dennis P. Townsend, and Erwin V. Zaretsky 

Lewis Research Center 

SUMAAARY 

A new rheological model incorporating four parameters: ambient viscosity, pres- 
sure viscosity coefficient, a lubricant factor, and transition shear stress was developed 
to describe the traction for a spinning point contact under elastohydrod 3 mamic conditions. 
The model is based on an exponential pressure viscosity model for shear stresses lower 
than 6.89x10 newtons per square meter (1000 psi). Above this shear stress, the expo- 
nential law is assumed to be followed until the ratio of shear /normal stress exceeds a 
certain value called the lubricant factor which is a pseudo -coefficient of friction. For 
all higher pressures, the shear stress is obtained by taking the product of the lubricant 
factor and normal stress. 

The analysis was compared with experimental spinning torque data reported in ref- 
erence 7 for a synthetic paraffinic oil, a di-2-ethylhexyl sebacate, a super -refined naph- 
thenic mineral oil, and a polyphenyl ether lubricant. Gtood correlation between theoret- 
ical and experimental results was obtained for these lubricants. The rheological model 
which limits the shear at high pressures but allows for low shear stresses in the micro - 
slip region, enables the spinning torque to be predicted to a sufficient degree of accuracy 
for most applications. 


INTRODUCTION 

In recent years, a complete evolution has occurred in ball bearing design and anal- 
ysis. Early work assumed that at the ball -race contact, spinning and rolling would occur 
at one raceway and only pvme rolling at the other raceway (refs. 1 and 2). This work wp.s 
based on the premise that there was a constant coefficient of friction at the ball -race 
contacts. Work reported in reference 3 and subsequently in references 4 and 5 indicated 
that the friction was not dependent upon a single value of friction coefficient but varied 



according to the magnitude of the contact stress, contact geometry, and lubricant type 
and viscosity. As a result, no simple analysis utilizing a constant value for the coeffi- 
cient of friction could be applied to a ball -race contact. Reference 6 recognizes this 
phenomenon and presents an analysis for the kinematics of a ball bearii^ considerir^ 
elastohydrodynamic effects and assumed lubricant rheological properties. However, 
these properties are not completely defined, resulting in some second -order inaccura- 
cies in the analysis. The authors of reference 7 developed a theoretical analysis based 
on the elastohydrodynamic theory of lubrication for spinning torque. This model pro- 
vided close agreement with the experimental data reported (ref. 7). 

In developing the analysis of reference 7 it was found that the assumption of a New- 
tonian fluid with an exponential pressure viscosity relation gave impossibly high values 
of torque. It was necessary therefore to introduce a cutoff point at which the pressure 
viscosity exponent decreased to a much lower value. The model of reference 7 yielded 
satisfactory results when the calculated shear rate was of the order of reciprocal 
seconds. Similar models were applied to the lubricants investigated experimentally in 
reference 3 and the results are reported in reference 8. 

Comparison with data of other researchers (ref. 9) showed the cutoff to occur at a 
much higher pressure when the film thickness was much greater than those of refer- 
ences 3 and 7. In other words, a lower shear rate resulted in a higher cutoff pressure. 

The objective of the work reported herein was to modify the composite viscosity 
model to take into account the shear rate dependency effects. The theory using the mod- 
ified model will be compared with the experimental results for a synthetic paraffinic lu- 
bricant reported in reference 7 and with the experimental results for a di-2-ethylhexyl 
sebacate, a super -refined naphthenic mineral oil, and a polyphenyl ether as reported in 
reference 3. 

The modified rheological model was suggested by remarks made by K. L. Johnson 
in a formal discussion presented in reference 7. He stated that ”at high pressvire, the 
film exhibits a critical shear stress which is approximately proportional to the pressure 
and decreases slightly with temperature , somewhat similar to a granular solid. This 
hypothesis leads to an approximately constant 'coefficient of friction’ t/s independent of 
film thickness, as would be observed in dry sliding or boundary friction. ” 


APPARATUS AND SPECIMENS 
Spinning Torque Apparatus 

4 

A spinning torque apparatus (see figs. 1(a) and (b)) described previously in refer- 
ences 3 and 4 was also used for the tests reported herein. The apparatus essentially 
consists of a turbine drive, a pneumatic load device, an upper and lower test specimen. 
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Figure L - Spinning-torque apparatus. 


3 



Front view 


r Contact ellipse 


\ 

' X 



Isometric view 

CO -9382 


(c) Lower test specimen. 
Rgure 1. - Concluded. 


a lower test -housing assembly incorporating a hydrostatic airbearing, and a torque- 
measuring system. In operation, the upper test specimen is pneumatically loaded 
against the lower test specimen through the drive shaft. As the drive shaft is rotated, 
the upper test specimen spins in the groove of the lower test specimen. This causes an 
angular deflection of the lower test -specimen housii^. This angular movement is sensed 
optically by the torque -measuring system and is converted into a torque value. During a 
test, the torque is continuously recorded on a strip chart. 


Specimens 

The upper test specimen is a conventional 12. 7-millimeter- (1/2-in. -) diameter 
bearing ball made of SAE 52100 steel. The ball has a nominal Rockwell C hardness of 61 

O 

and a surface finish of 5x10" meters (2 pin. ) rms. The lower test specimen (fig. 1(c)) 
is a 12. 7 -millimeter - (1/2 -in. -) diameter ball from the same heat of material as the up- 
per test specimen. The lower ball is modified by grinding a flat on one side and a cylin- 
drical groove of radius Rq (fig. 1(c)) on the other. The groove simulates the race 
groove of a bearing. The axis of the groove is parallel to the flat. The groove radius 
expressed as a percentage of the upper -ball diameter is defined as the ball -race con- 
■ formity. The specimens used in these tests were groimd to ball -race conformities of 51, 
55, and 60 percent. The simface finish of the cylindrical groove was approximately 
4 5X10"® to 15x10"® meter (2 to 6 pin.) rms. 
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THEORETICAL ANALYSIS 


An elastohydrodynamic analysis for the determination of the torque of a ball spinning 
in a nonconforming groove, shown in figure 1, was developed in reference 7. In the anal- 
ysis of reference 7, the area of contact is elliptical, as shown in figure 2(a). In this 
ellipse an elastohydrodynamic film can theoretically be predicted, due to the spinning ve- 
locity, except within the inscribed circle of radius b. The frictional force over the 
ellipse is determined by making the system into a number of elemental rollers of width 
dy as shown in figure 2(b), and determining the frictional force on each roller. These 
rollers slide into and out of the contact region as the major axis of the ellipse is crossed. 

The torque due to the region outside the inscribed circle and inside the contact 
ellipse is computed by integrating the moments of the elemental rollers over this region; 
that is. 



( 1 ) 


To determine the frictional force on each elemental roller, it is split into elements 
of width dx. The friction force on each elemental roller is thus given by 


dF = 



( 2 ) 


where t is shear stress. For a Newtonian fluid, assuming a linear velocity gradient 
and ignoring the velocity in the x direction which is small, the approximate shear stress 
is given by 


(3) 

h 

where h, the film thickness, is given by reference 10 for a roller and modified for each 
elemental roller by 


h = 


1. 6q!°' 


,0.7 


W 


0. 13 


(4) 
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(a) Contact ellipse for ball in nonconforming groove. 


— -j |— dy 



(b) Elemental roller. 


Figure 2. - Ground lower ball showing contact ellipse and elemental roller. 
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where 


For steel, E* is 22.3x10^® newtons per square meter (32. 3x10® psi). 

The assumption of a Hertzian stress distribution leads to a load distribution W; that 
is, 



( 6 ) 


where P is the normal load. 

Using the exponential pressure viscosity relation, the viscosity at any location x, y 
would be given by 


fi = 



(7) 


where S is the contact pressure at that point. If a Hertzian distribution is assumed, the 
contact pressure is given by 


1. 5P 
7rab 



( 8 ) 


The proposed rheological model is based on the assumption that the foregoii^ rela- 
tions are valid at low shear stresses but a point is reached above which the shear stress 
becomes proportional to the normal stress. These conditions may be represented by the 
following relations: 


T = 


0'S (oy 


(9a) 


T = fS 


and -^ > fS 

° h 


(9b) 
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The introduction of a lubricant factor f serves to limit the shear stress at high 
pressures and shear rates to a fraction of the normal stress. However, there are cases 
at low pressures where the shear stress is greater than fS. Where this occurs, it is 
more reasonable to use the viscous model. Therefore, a transitional shear stress 
can be introduced to allow for large values of shear stress at low pressure where the 
viscous model shows a shear stress greater than fS. Thus, an additional condition is 
imposed wherein equation (9a) should be used where 

T < T and T < T < fS 
c c 

In order to visualize the behavior of the fluid under the foregoing conditions, figure 3 
shows a plot of shear stress as a fvmction of normal stress (pressure) for a synthetic 
paraffinic Ivibricant at various shear rates. For purposes of comparison, the shear 
stresses computed by equations (9a) and (9b) may be substituted into equation (3), and an 

50x10* 



Pressure (normal stress), psi 


Figure 3. - Shear stress plotted against pressure for 
synthetic paraffinic lubricant at various shear 
rates. 


8 



equivalent viscosity may be obtained for any shear rate. Under the method of reference 7 
this form of the equivalent viscosity would be used to calcxdate the spinning torque. Fig- 
ure 4 shows this equivalent viscosity for various shear rates and also a comparison with 
the previous rheological model from reference 7. The previous model is seen to be ap- 
proximately equivalent to a special case of the newer, more general, model. 

Computation of the moment due to the outer region of the ellipse is undertaken by us- 
ing equations (1) and (2). Equation (4) is used to compute the film thickness. The rhe- 
ological relations (eqs. (9a) and 9(b)) are used to compute the shear stress at each point. 
By combining and integrating these equations twice numerically, the moment about the z 
axis is obtained. 

There remains the area within the inscribed circle of radius b where a conventional 
elastohydrodynamic film theoretically is impossible to maintain. However, separation of 



Pressure, psi 

Figure 4. - Equivalent viscosity plotted against pressure using new rheological 
model for synthetic paraffinic lubricant at various shear rates. 
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the surfaces is made possible by means of a microasperity elastohydrodynamic film 
(ref. 11) or squeeze film. Because the area outside the inscribed circle from y equals 
-b to +b is small, it is neglected. If the inscribed circle is then divided into elemental 
rings of radius r and width dr, the torque Mg for the inscribed circle becomes 


Mg = 2t7 




r^ dr 


The moment over the inscribed circle is then 


Mg = 


27rT^b‘ 


(10) 


( 11 ) 


In addition to the effect of the lubricant within the Hertzian contact region, the effect 
of the viscous drag of the lubricant outside this region must also be considered (fig. 5). 
The expression for this moment is given in reference 5 as 
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This may be integrated numerically over the region outside the contact ellipse. The total 
spinning torque is therefore the sum of the moments given by equations (1), (11), 
and (12) 


Ms = + Mg + Mg 


(13) 


RESULTS AND DISCUSSION 


The spinning moments for a ball spinning in a nonconforming groove with a synthetic 
paraffinic lubricant and conformities of 51, 55, and 60 percent were computed using the 
revised rheological model and compared with the experimental values from reference 7. 
The theoretical results which are given in figure 6 show good agreement with the exper- 
imental data of reference 7 which are shown for comparative purposes. The lubricant 
parameters used in the rheological model were as follows: 


Ambient viscosity, N-sec/m^ (Ib-sec/in. ^) 0.414 (6x10”^) 

Pressvire viscosity coefficient, a, m^N (psi”^) 1. 33x10"® (0. 92x10"^) 

Transitional shear stress, t^, N/m^ (psi) 6.89X10® (1000) 

Lubricant factor , f 0.07 


(The preceding values for viscosity and pressure viscosity coefficient were the same as 
reported in ref. 7. ) 

For the di-2-ethylhexyl sebacate Ixibricant tests which were reported in reference 3, 
the following fluid parameters were used: 


2 2 

Ambient viscosity, u , N-sec/m (Ib-sec/in. ) 

^ 2-1 
Pressure viscosity coefficient, a, m /N (psi” ) 

9 

Transitional shear stress, t^, N/m (psi) . . 
Lubricant factor, f 


. 0.016(0.24x10."®) 

1.45X10"® (l.OxlO"^) 
. . 6. 89X10® (1000) 

0.045 
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40 60 80 100 120 140 160 180 200220x10^ 

Maximum Hertz stress, psi 

Figure 6. - Comparison of experimental and theoretical 
spinning torque for synthetic paraffinic lubricant for 
conformities of 51. 55, and 60 percent. Ambient 
viscosity Mg = 0.414 N-sec/m^ (6xlCr Ib-sec/in. 
pressure viscosity coefficient a= 1. 33x10'° m^/N 
(0.92x10'^ psi'*); transitional shear stress Tg = 
6.89x10® N/m^f 1000 psi); lubricant factor f=0.07; 
spinning speed. 1050 rpm. 
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The value of a for the di -2 -ethylhexyl sebacate is that reported in reference 12. The 
experimental and theoretical cvirves which are shown in figure 7 are in good agreement. 

For the super -refined naphthenic mineral oil, good agreement was obtained between 
the experimental results from reference 3 and theoretical results using the rheological 
model with the following parameters: 


Ambient viscosity, N-sec/m^ (Ib-sec/in. 0.079 (1.14x10“^) 

Pressure viscosity coefficient, a, m^/N (psi“^) 2. 15x10"® (1. 48x10“^) 

Transitional shear stress, t^, N/m^ (psi) 6.89x10® (1000) 

Lubricant factor, f . . . 0.075 


The value of a used here is the same as reported in reference 13. These results are 
shown in figure 8. 



Maximum Hertz stress, psi 


Figure 7. - Comparison of experimental and theo- 
retical spinning torque for di -2 -ethyl hexyl 
sebacate. Ambient viscosity Pp = 0.016 N -sec/ m‘ 
(0. 24x10'^ Ib-sec/in. \ pressure visrosity co- 
efficient a= 1.45x10”° m‘/N (1.0x10'^ osi'ljj 
transitional shear stress Tj. = 6.89x10° N/m^ 
(1000 psi); lubricant factor f = 0.045; spinning 
speed, 950 rpm; conformity, 51 percent. 
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Maximum Hertz stress, psi 

Figure 8. - Comparison of experimental and theo- 
retical spinning torques for super -refined 
naphthenic mineral oil. Ambient viscosity \ 1 q = 
0. 079 N-sec/m^ (1. 14x10’^ Ib-sec/in. ‘); pres- 
sure viscosity ^efficient a = 2. ISxlO'* m’/N 
(1.48xl0‘‘* psi"*); transitional shear stress 
= 6. 89x10^ N/m^ (1000 psi); lubricant factor 
f= 0.075; spinning speed, 950 rpm; conformity, 
51 percent. 


The comparison between the experimental (ref. 3) and theoretical results for the 
polyphenyl ether is shown in figure 9. While there was reasonably good correlation be- 
tween the theoretical and the experimental results, it was not possible to obtain as good 
a correlation as with the other fluids even though a number of different lubricant param- 
eters were tried. The parameters giving the best correlations are as follows: 

2 2 S 

Ambient viscosity, u, N- sec/m (Ib-sec/in. ) 0.8(11.6x10") 

^21 ft 4 

Pressure viscosity coefficient, a, m /N (psi" ) 4.64x10" (3. 2x10" ) 

Transitional shear stress, N/m^ (psi) 6.89x10^(1000) 

Lubricant factor, f 0. 07 
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Maximum Hertz stress, N/m^ 


I I \ I I I I L_J , 

60 70 80 90 100 110 120 130 140x10^ 

Maximum Hertz stress, psi 

Figure 9. - Comparison of theoretical and experi- 
mentai spinning torques for polyphenyl ether 
(5P4EI. ^bient visOTSity Pu" 0.8N-sec/m‘ 
(11. 6x10'^ Ib-sec/in. ^)- pressure viscosity 
coefficient 0=4.64x10"® mWN 13.2x10''’ psi ‘H; 
transitionai shear stress T(. = 6. 89x10° N/m‘ 
(1000 psi); lubricant factor f=0.07; spinning 
speed, 950 rpm; conformity, 51 percent. 


(The preceding value of a is that given in ref. 14. ) 

CONCLUDING REMARKS 

The results of this study have shown that the traction in a sliding elastohydrodynamic 
contact may be computed on the basis of four fluid parameters. Of these parameters 
only the following three are critical: 

(1) Ambient viscosity 

(2) Pressure viscosity coefficient a 

(3) Lubricant factor f 

The fourth parameter, the transitional shear stress t , is merely a device to allow the 

In' 
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Newtonian viscosity relation to be applied at low pressure when the ratio t/s is greater 

than f and may also be greater than unity. An identical value for the transition shear 
fi 2 

stress (6. 89x10 N/m (1000 psi)) was used for all lubricants in this study. 

The pressvme viscosity coefficient used in the analysis was, in each case, taken 
from published values. The lubricant factor was adjusted imtil good correlation was ob- 
tained. In general, a lubricant factor in the range 0. 07 to 0. 075 gave good correlation. 
The exception to this was the di-2-ethyhexyl sebacate for which a lubricant factor of 

0. 045 was found. The rheological model enables the spinning torque to be predicted to a 
sufficient degree of accuracy for most applications. 


SUMMARY OF RESULTS 

A new rheological model incorporating four parameters: ambient viscosity, 
pressure -viscosity coefficient, a lubricant factor and transition shear stress was devel- 
oped to describe the traction for a spinning point contact under elastohydrodynamic 
conditions. 

The new rheological model is based on an exponential pressure -viscosity model for 
shear stresses lower than 6. 89 newtons per square meter (1000 psi). Above this shear 
stress, the exponential law is assumed to be followed until the ratio of shear /normal 
stress exceeds a certain value called the lubricant factor which is a pseudo -coefficient of 
friction. For all higher pressures, the shear stress is obtained by taking the product of 
the lubricant factor and normal stress. 

The analysis was compared with experimental spinning torque data reported in ref- 
erences 3 and 7 for a synthetic paraffinic oil, a di-2-ethylhexyl sebacate, a super -refined 
naphthenic mineral oil, and a polyphenyl ether lubricant. The following results were 
obtained: 

1, Good correlation between theoretical and experimental results was obtained for 
the synthetic paraffinic oil, the di-2-ethylhexyl sebacate, the super -refined naphthenic 
mineral oil, and the polyphenyl ether lubricant. 

2. The rheological model which limits the shear at high pressures but allows for low 
shear stresses in the microslip region, enables the spinning torque to be predicted to a 
sufficient degree of accuracy for most applications. 

I 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, February 20, 1973, 

501-24. 
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APPENDIX - SYMBOLS 


a major semiaxis of contact ellipse, m (in.) 

b minor semiaxis of contact ellipse, m (in. ) 

b’ semiwidth of contact ellipse at y, m (in.) 

E modulus of elasticity, N/m (psi) 

E’ materials properties factor, N/m^ (psi) 

F friction force 

f lubricant factor 

h film thickness, m (in. ) 

h|^ value of h at x = b, m (in. ) 

h^ minimum distance between ball and groove, m (in. ) 

K constant definir^ outer bovindary of integration 

M total spinning torque in Hertzian ellipse, N-m (lb -in. ) 

Mj spinning torque in Hertzian ellipse outside inscribed circle, N-m (lb -in. ) 

Mg spinning torque in Hertzian ellipse inside inscribed circle, N-m (lb -in. ) 

Mg spinning torque due to viscous drag outside Hertzian ellipse, N-m (lb -in. ) 
P normal load, N (lb) 

R radius of equivalent cylinder, m (in.) 

Rg radius of ball, m (in. ) 

Rq radius of groove, m (in. ) 

r polar coordinate 

Tq value of r at outer boimdary of Hertzian ellipse, m (in. ) 

S contact stress, N/m (psi) 

u relative velocity, m/sec (in. /sec) 

W load per unit width, N/m (Ib/in. ) 

x,y,z Cartesian coordinates, m(in. ) 
a pressure viscosity exponent 

2 2 

p absolute viscosity, N-sec/m (lb -sec/in. ) 

2 2 

Pq ambient viscosity, N-sec/m (lb -sec/in. ) 
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a Poisson's ratio 

2 

T shear stress, N/m (psi) 

2 

value of T at x = b, N/m (psi) 

2 

Tp transition shear stress, N/m (psi) 
^ polar coordinate 
CO angular velocity, rad/sec 
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CONTRACTOR REPORTS: Scientific and 
technical information generated under a NASA 
contract or grant and considered an important 
contribution to existing knowledge. 


TECHNICAL TRANSLATIONS: Information 
published in a foreign language considered 
to merit NASA distribution in English. 

SPECIAL PUBLICATIONS: Information 
derived from or of value to NASA activities. 
Publications include final reports of major 
projects, monographs, data compilations, 
handbooks, sourcebooks, and special 
bibliographies. 

TECHNOLOGY UTILIZATION 
PUBLICATIONS: Information on technology 
used by NASA that may be of particular 
interest in commercial and other non-aerospace 
applications. Publications include Tech Briefs, 
Technology Utilization Reports and 
Technology Surveys. 

publications may be obtained from: 


Details on the availability of these 

SCIENTIFIC AND TECHNICAL INFORMATION OFFICE 
NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C. 20546 



